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Abstract 

We study a lattice sigma model which is expected to reflect the 
Anderson localization and derealization transition for real symmetric 
band matrices in 3D. In this statistical mechanics model, the field 
takes values in a supermanifold based on the hyperbolic plane. The 
existence of a diffusive phase in 3 dimensions was proved in [2] for 
low temperatures. Here we prove localization at high temperatures 
for any dimension d > 1. Our analysis uses Ward identities coming 
from internal supersymmetry. 



1 Introduction 

It is well known that the study of localization properties in a disordered 
material can be translated to the study of correlation functions in a lattice 
field theory, with an internal hyperbolic supersymmetry (SUSY), [9j El El [8]. 
In the physics literature one usually assumes the sigma model approximation, 
which is believed to capture the essential features of the energy correlations 
and transport properties of the underlying quantum system. 
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The SUSY field theories which are equivalent to the Anderson tight- 
binding model and random band matrices are difficult to analyze with math- 
ematical rigor in more than one dimension. In this context, Zirnbauer intro- 
duced a lattice field model which may be thought of as a simplified version of 
one of Efetov's nonlinear sigma models [TO], [3]. In Zirnbauer's sigma model 
the field takes values in a target space if( 2 l 2 ) which is a supermanifold exten- 
sion of the hyperbolic plane. The model is expected to reflect the spectral 
properties of random band matrices, such as localization and diffusion, in 
any dimension. In [10] localization was established in one dimensional chain 
by analyzing the transfer matrix. We refer to [2] for a historical introduction 
and motivations. 

More recently the existence of a 'diffusive' phase at low temperatures {(3 
large) has been proved for the if( 2 l 2 ) model in three or more dimensions, 
see [2]. For (3 small, a localized phase was expected. However, unlike con- 
ventional statistical mechanics models, the if( 2 l 2 ) model has a noncompact 
hyperbolic symmetry and so high temperature expansions cannot be done 
in the usual way. In fact, it is known that the bosonic hyperbolic sigma 
model in 3D has no localized phase because its effective action is convex for 
all (3 > 0. On the other hand, numerical simulations [I] indicated that the 
SUSY hyperbolic sigma model has a phase transition for (3 < (3 C ~ 0.038. 

In this paper we show that for any dimension d > 1 the if( 2 l 2 ) model 
exhibits localization for (3 l / 2 \n(3~ l < l/(2d— 1). Thus the sigma model 
approximation captures the physics of both localization and diffusion. More- 
over, for a one dimensional chain we recover localization for all values of (3. 
Localization is also expected in 2D (see [2] section 1.4 and 4.3), for all val- 
ues of (3 by both the renormalization group and by a simple saddle analysis. 
However, a rigorous proof is still missing for this case. 

The techniques employed in this work to prove localization are quite dif- 
ferent from the ones used in [2] to prove extended states. The two papers can 
be read independently. The only common point is the use of supersymmetry 
to prove some identities. In the present case supersymmetry is applied only 
to prove that the partition function is normalized to 1. We refer to Sec- 
tion 4 and Appendix C in [2] for an introduction to supersymmetric Ward 
identities. 
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1.1 The model 



Let A be a finite subset in Z d and tj a real variable for each site j G A. We 
will consider periodic or Neumann boundary conditions. Accordingly for any 
two points x, y G A, \x — y\ will denote the Euclidian distance on the lattice 
(Neumann be) or the periodized lattice (periodic be). 
We introduce the probability measure 



dff A {t)= H^e- F ^) e -MUt) x 



det 



where dtj is the Lebesgue measure, F is the kinetic part and M is the mass 
term: 

F A (Vt) =/3 ^ (coshfo - fy) - 1) 

Ml(t) = J^Ejicoshtj - 1) . (1.2) 

jeA 

We denoted by (jf) the nearest neighbor pairs \ j—j'\ = 1. F) £ K {t) is a positive 
definite matrix defined by 

M)« =o K-j|>i 

{pftii =-(3 \i-j\ = l (1-3) 

(DDjj =P[2d+V J )+e j e-*' i = j, 

v i = E [ e ' fc " 4j - !] ( L4 ) 

fc,(jfc) 

and £j > are regularizing parameters that are necessary to make the integral 
well defined. We remark that V 3 ■ + 2d > for all t configuration and all j. 
Finally (3 > is a parameter that can be interpreted as a measure of the 
temperature or the disorder. 

Note that the definition of the matrix D differs from the one introduced 
in [2] eq. (1.1). It is not difficult to see that when you mix the term ^2 k tk 
and the determinant in the effective action (1.2) of [2] the result is indeed 
the determinant of the matrix D above, more precisely 

|i-j|>l 
Aij = e ti D ij e ti = { -(3e u+t > ji-jj = 1 (1.5) 

£*,(*') e*"* *=3 
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where A is the matrix introduced in [2] eq. (1.1). For technical reasons the 
above representation is more convenient when we want to prove diffusion as in 
[2] while the other is more practical when we study localization (except in the 
proof of Theorem 2 where we will go back to the "diffusion" representation 
for a while). 

Now, for any function f(t) we will define its average by 

(/(*)> = / dfi £ A (t) f(t) . (1.6) 

Normalization and choice of Ej By internal supersymmetry (see [2] Sect. 
4 and eq.(5.1)) this measure is already normalized to 1 so the partition func- 
tion is 

Z% = J dff A (t) = 1 • (1.7) 

This identity is true regardless of the boundary conditions and the values of 
(3 or Sj as long as the integral is well defined. Since we consider (3 > fixed 
we only need Sj to be non zero at one lattice point. In the following we will 
consider three cases. 

1. Uniform pinning: Ej = e < ± for all j G A. The measure is translation 
invariant with periodic be. The correlation function in this case has a 
divergent prefactor l/e in the localized regime. 

2. Two pinnings: e x = e y = 0(1) and Ej = for all other points. This is 
the analog of inserting two electrical contacts in a metal sample. 

3. One pinning at j = 0: 6q = 0(1) and Ej = for all other points. This is 
more suitable for an interpretation of the model as a random walk in a 
random environment. Our results suggest that edge reinforced random 
walk (see [7]) will also localize when the reinforcement is strong. 

The observable. We will study the correlation function 

G xy := /;,,; (1.8) 

where x, y can be any two points on the lattice such that both e x > and 
E y > 0. This observable does not give information on localization properties 
in the case of one pinning point. In such case a good observable to study is 

Oj = e +t > /2 (1.9) 
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where j is any point in the lattice. This observable is analogous to xj^ in 
the notation of [7] where e is an an edge (j,f)- 

1.2 Main results 

Theorem 1 Let e x > 0, e y > and X^eA £ j < 1- Then for alio < (3 < (3 C 
(P c defined below) the correlation function G xy (11.81) decays exponentially with 
the distance \x — y\. More precisely: 

(G xy ) < Coie^ + e- 1 ) [l,e^ c d ] {x ~ vl , (1.10) 
where Cd = 2d — 1 , Cq is a constant and 

— f 00 rlf 



. 27T 

Finally (3 C is defined by: 

Ip c e^-% d =l . (1.12) 
Our estimates hold uniformly in the volume. 

Remark 1 The integral Ip satisfies Ip < 1 V/3 > 0: 

( /*oo ( /*oo 

Ip=sJ^ / cfte"^ 0811 ^ < \jt dt cosh(t/2) e -0(«***-i) = 1. 

(1.13) 

More precisely 

r /j < o.i5 

^ ~ \ ce~^ /3 » 1 V ; 

where c > 1 is a constant. 

Remark 2 The constraints e x > and > exclude the case of one 
pinning. Moreover YljeA £ i — ^ i m phes e < m when = e is constant. The 
case of one pinning is covered by Theorem 2 below. 
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Main consequence For d = 1 the critical beta is (3 C = oo since 

IpePte-Va = Ip < 1 V/?>0. (1.15) 

Therefore the correlation function decays exponentially for all values of f3. 
On the other hand for d > 1 we obtain localization only for small j3 since 
& < (2d- I)" 2 < 1. 

Theorem 2 Let e = 0(1) and = Vj ^ 0. T/ien /or all < (3 < (3 C 

(P c defined below), the field t x wants to be as negative as —\x\. More precisely 
(O x ) (defined in (11.91) ) decays exponentially with the distance \x — y\ 

(O x ) < Co [lpe^-V Cd ] ]x] , (1.16) 

where q and Ip are defined in Theorem 1 above and Co is a constant. Finally 
(3 C is defined by: 

Our estimates hold uniformly in the volume. 

Main consequence. For d = 1 the critical beta is (3 C = oo since 

Ipe^ Ci ~ x) c d = Ip<\ W(3>0. (1.18) 

Therefore (O x ) decays exponentially for all values of f3. On the other hand 
for d > 1 the result holds only for small (3 since (3 C < (2d — 1)~ 2 < 1. 

Acknowledgments. It is our pleasure to thank A. Abdesselam for discus- 
sions and suggestions related to this paper. A special thanks to M. Zirnbauer 
who explained the model to us and shared his many insights. 

2 Proof of Theorem H 

We want to estimate 
(G xy ) = J d^{t)G xy = J J] jLe-^) e -^) x ^ detD e A{t) D -i . 



(2.1) 

The proof in done in four steps. 
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Step 1. We mix the observable G xy and a piece of the probability measure 
namely V det D. The key identity is 



Met D\{t) D~y = y/D£ \j D~y det D%(t) (2.2) 
(remember that D~y > 0). The first term is bounded by 

This is proved in Lemma [H Inserting this in (12 .ip we have 

r i / e t x /2 e t y /2\ 

(G xy ) < J dS A (t)y/[D£d4tD%{t)]^— + —) (2.4) 



where 

jeA v 27r 

Unlike the measure <i/i^ given by (11. II) . this measure is no longer normalized 
to 1. 

Step 2. We need to extract some decay in the distance \x — y\. This is 
hidden in D~ y l det-D. By some combinatorial arguments (the proof is given 
in Lemma [2] ) we can write 

[D^detDUt)} = J2? hl detD U ( 2 - 6 ) 

where the sum ranges over non self intersecting paths 7 made of nearest 
neighbor pairs in A starting at x and ending at y. Let I7I, denote the length 
of 7 and let A 7 be the corresponding set of lattice points and set = A\A 7 . 
Finally D\ c is the matrix one obtains by deleting the rows and columns 
corresponding to the lattice points j G A 7 . It is is exactly like the matrix 
D £ Al but defined on the complement of 7, A 7 and with modified masses: 

Dfj =0 |i-j|>l 

D !j =~P l*-il = l (2.7) 



2d + V l 



+ eie tl % = j 
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where i,j G A 7 and 



V i= [e tk ~ h -l] (2.8) 

fceA 7 ,(fc«) 

Si = Si + (3 e<fc • ( 2 - 9 ) 

fceA 7 ,(fci) 



By combining (12.41) and (12.61) we have 

(G xy ) < Ys^' 2 J (7= + ^) \/^D|7 • (2.10) 

Step 3. The measure di>^(t) defined in (12.51) can be factored as a measure 
on A 7 times a measure on the complement set A 7 

dul(t) = dul^t) du^(t) e~ F ^t) ; ( 2.ii) 

where 

^ 7 (V£) = Yl (3(cosh(t 3 -t k )-l) (2.12) 

(j,fc),fceA 7 j^A 7 

describes the interaction between A 7 and A°. Then the integral in (I2.10p can 
be written as 

/• / P W2 P *j//2\ , r / e t*ri e t y /2\ 

.1 d " m ( + v^) = y ^ w ( + v^) z ' k(h) ' 

(2.13) 

where we defined 

ZZcfo) = | ^detZ^ e"^^ 

= / n 72= e_ ^ (vt) e_M ^ (t) e_Fs7(Vi) • (2 - i4) 



j6A 7 

Note that Z\ c (i 7 ) is still a function of the t variables along the path {tfc}fceA 7 
(they are not integrated). Now Z\ c (t^) is almost equal to the partition 
function 



i=z u= J ^y*) = / du u^\i deiD u 



II 4= e^ (V V M ^ W ^/det^le , (2.15) 

jeA 7 v 
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where Z Ac = 1 by supersymmetry (see (11.71) ). Comparing (12.141) and (12.151) 
we see that there are two main differences: 

- the mass term M Ac in (12.141) depends on e instead of e and is smaller 
than what it should be 

M E A c < M £ A c 

since M £ contains additional mass terms; 

- the exponent in (I2.14p contains the additional factor — Fg 7 (Vt) coming 
from the kinetic interaction between points on A 7 and points on A T . 

This last term is helping us since it makes the integral smaller. We will use 
it to recover the missing mass. This is done in Lemma [3] below. The key 
ingredient is a global translation on the t variables. The result is 

(G xy ) < e 1 £/jH/VI»rl / d^{t) + (2 .i 6) 

where \d^\ < {2d— 2)\^\ +2 is the number of points inside A 7 on the boundary 
with A 7 . 

Step 4. We are left with an integral along the path 7. The integral in 
( I2.16P is bounded by 

/ e W2 e ty/2\ f jx jy \ 

dvut) e — + e —) < 4= + ^= 4 7 ', 



where 



/; = / _iL e */2 e - £ ,(cosh 4 -i) = J_ (217) 

-00 \J2tx \f£^ 



00 



dt «/„„„u + ii 1 



h = I —e-^ ht -V = ^j- 2 I, (2.18) 



oo 
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where Ip was defined in ( I2.17p . In the same way If = l/^/Sy~. Inserting all 
this we have 

(G xy ) < jfL + LWefiW^p (2.19) 

^ X n>\ X -y\ 

<Co(- + -) [e^I,c d f- yl 



where c d = (2d— 1), Co is a constant and the second inequality holds since the 
number of self-avoiding walks made of n steps is bounded by 2d(2d— l) n < 2c^ 
and \dj\ < (2d — 2)\ / y\ + 2. Finally the sum over n is convergent since 
(e^ 2d -Vl(3C d ) < 1. ' 

This concludes the proof of Theorem 1. □ 

2.1 The lemmas 

Lemma 1 The following inequality holds: 

ptx pty 

D-y 1 < — + — • (2.20) 



&x &y 



Proof By Cauchy-Schwartz inequality 



Since (/, Df) > V. Eje~ tj f 2 for any fj eKwe have 

e x e tx 

Hence the result. □ 
Lemma 2 For any invertible matrix M on A we have the following identity: 
[M-MetM] = \{-M xn ){-M nn )...{-M ]my )] det A ,M (2.21) 

Jx V =(jl,—jm.) 

where 7 is any non self intersecting path starting at x and ending at y. 
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Proof This is a classical formula arising from the fact that every permu- 
tation can be decomposed as a product of cycles. One may derive it using 
the representation of a determinant as a sum over a gas of disjoint non self 
intersecting closed paths: 



detM= A(Lx)---A{L % 



L\,...Lp 

where a loop L = (jx, . . . , j m ) is an ordered set of m distinct points and 
A(L) f - [(-M jlj2 )(-M j2h ) ■ ■ ■ (-M jmjl )] m > 1 

The sign (— l) m_1 is the number of pairs inside the loop that need to be ex- 
changed in order to recover the trivial permutation. Now since [M^detM] = 
»d detM, the derivation selects only loops that contain the pair yx. The 

ONI y x 

corresponding matrix element disappears and the loop becomes a path from 
x to y. The sign —1 from —M yx cancels the global —1 in front of the product. 
□ 

Lemma 3 For any configuration of {t/~\ k e A 7 } ; the conditioned partition 
function Z\ c (t 7 ) given by (I2.14j) is bounded by 

Zl.it,) < e ^ 6 A^ 1 - etfc - t V^^ 1 ~ e "* ) < e^V^I , (2.23) 
where t* is any real number satisfying 

t* > , and t* >t k \/k G A 7 , (2.24) 
and dl is the number of points nearest neighbor to k that do not belong to 

^ = #{j£A 7 | \j-k\ = l}. 



A 7 ; 



Proof Before doing any bound we perform a global translation inside the 
integral: 

tj -> tj + 1* Vj G A 7 . (2.25) 
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Then inside the exponential we have: 
F A c(Vi) i- F A c(Vt) 

F dj (Vt) i-> F dj (Vt + f)= Yl $ (cosh^ +t* - t k ) - r 

(i,A;),fceA 7 ,ieA= 



(2.26) 
(2.27) 

(2.28) 



In the numerator we have 



det D £ A 



det DlC 



so the only effect of the translation is to modify the mass £j defined in (12.91) 
at each lattice point. After the translation 



7 7 



n 



dt , 



e-^ (Vt \/detmr*e- M ^ ( * +r) e-^( v ^) 



^2.29) 



dt , 



II \ : 



e -F A c(V*) e -M|| (*) ^ detD ee-** e M % (t)-M^(t+t*)-F d ,(Vt+n 



where we used Z £ £ ' = 1 (see (11.71) ) and we defined 



^2.30) 



£ rrj = £ ie "~**(cosht,-l) - ej(cosh(ij+t*)-l) - Y (3 (cosh(t j+t* -t k )-l). 

fceA 7 ,(j,fc) 

(2.31) 

To conclude we shall prove that the right hand side of f!2.29|) is bounded by 
the rhs of (122511 : 

i£A 7 

We distinguish two cases. 
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Case 1. When j is far from the path 7 that is \j — k\ > 1 for all k G A 7 
then £j = Ej and we have 

£rrj = Ej [e"**(cosh^- - 1) - (cosh(t j + £*) - 1)] 

= Ej [e f j sinh(-f) + 1 - e - **] < Ej(l - e~ 1 *) , (2.32) 

where the last inequality holds for t* > 0. 



Case 2. When j is near to the path 7 that is \j — k\ — 1 for some fce A 7 
then ^ = £j + Efc 6 A 7 ,(j,fc) ^ and we have 

£ rrj = £j [e^^coshtj - 1) - (cosh(t j + £*) - 1)] (2.33) 

+ <^E [ e * fe_< * ( cos H' - !) - (cosh(* 3 - + t* - t k ) - 1)] 

fceA 7 ,(j,fc) 

= £j - [e^ sinh(-r) + 1 - e~**] + [e^ sinh(£ fe - f ) + 1 - e^] 

<E,(l-e- tt )+(3j2 (1-e^*), 

fceA 7 ,(j,fc) 

where the last inequality holds if t* > and t* > VA; e A 7 . 
Finally 

J>rr,< [e J (l-e-n+ P £ (1 - e^) ] 

jeA 7 jeA= fceA 7 ,(j,fc) 

</?E ^ - ^ + E - e_r ) • ( 2 - 34 ) 
This concludes the proof of the lemma. □ 



3 Proof of Theorem 2 

The proof of Theorem 2 is almost identical to that of Theorem 1. This time 
there is no term D~y ensuring we can extract a path 7 connecting x to y. 
On the other hand, since there is a pinning only at one position the matrix- 
tree theorem (see p] for a simple proof and many references) applied to the 
"diffusion" representation A given in (I1.5P of the matrix D gives 

deU = Eoe^ II MaO = E II [^ +V ] > 
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where the sum is over the spanning trees on A made of nearest neighbor pairs 
(since Ajf = when \j — j'\ > 1). Therefore each term in the sum contains 
a path 7 from to x. Actually using (13. 1 p it is easy to see that 

A,} detA = — — detA (3.2) 

u e e to 

for all x G A, x ^ 0. Therefore Boe^A^ = 1 and 
det£> = e e t0 A^ detD = e e~ tx D^deW = e e~ tx ^ p hl det A oD , (3.3) 

70a; 

where Aq x = e~ to D^e~ tx and in the last term we applied Lemma[2j Inserting 
this result in (11.11) we have 



(O x ) = j dfx%(t) O x = J dul(t)^tDle^ 2 



(3.4) 



^A2 



= v^) / du £ A (t)e u/2 e- tx/2 Ij^f 3 ^ detD 

J V 70a; 

70a; 70a; 

< J2 ^ ^ $ < Co (c d e^I^ , 



7().r 



where dv £ K {t) was defined in (12T5I) . ZjLfa) in fl2TT4j) . ^ in 051) and the 

same definition holds for I £ . In the second line we used detA^-D = detD A c 

(see (12.71) ). We used Lemma 3 eq. (12.231) to bound Z]^ c {t y ). Finally the last 

inequality holds since (Qe^ -1 /^) < 1. This concludes the proof of Theorem 2. 
□ 
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